RANDOM QUANTUM CHANNELS II: ENTANGLEMENT OF RANDOM 
SUBSPACES, RENYI ENTROPY ESTIMATES AND ADDITIVITY 

PROBLEMS 



BENOIT COLLINS AND ION NECHITA 

Abstract. In this paper we obtain new bounds for the minimum output entropies of ran- 
dom quantum channels. These bounds rely on random matrix techniques arising from free 
probability theory. We then revisit the counterexamples developed by Hayden and Winter to 
get violations of the additivity equalities for minimum output Renyi entropies. We show that 
random channels obtained by randomly coupling the input to a qubit violate the additivity 
of the p-Renyi entropy, for all p > 1. For some sequences of random quantum channels, we 
compute almost surely the limit of their Schatten Si Sp norms. 



1. Introduction 

The relationship between random matrix theory and free probabihty theory Ues in the as- 
ymptotic freeness of random matrices. Asymptotic freeness, as it was discovered by Voiculescu 
(see e.g. [23]), usually predicts the asymptotic pointwise behavior of joint non-commutative 
moments. However in some cases it can also predict more. For example, in the case of i.i.d. 
GUE random matrices, it was showed by Haagerup and Thorbj0rnsen that even the norms 
have an almost sure behavior predicted by free probability theory. 

Quantum information theory is the analogue of classical information theory, where classical 
communication protocols are replaced by quantum information protocols, known as quantum 
channels. Despite the apparent simplicity of some mathematical question related to informa- 
tion theory, their resistance to various attempts to (dis)prove them have led to the study of 
their statistical properties. 

The Holevo conjecture is arguably the most important conjecture in quantum information 
theory, and the theory of random matrices has been used here with success by Hayden and 
Winter [TU [16] to produce counterexamples to the additivity conjecture of Renyi entropy for 
p > 1. These counterexamples are of great theoretical importance, as they depict the likely 
behavior of a random channel. In [12j, Hastings gave a counterexample for the case p = 1. It 
is also of probabilistic nature, but uses a very different and less canonical measure. 

In our previous paper [8], we introduced a graphical model that allowed us to understand 
more systematically the computation of expectation and covariance of random channels and 
their powers. In particular we studied at length the output of the Bell state under a random 
conjugate bi-channel and obtained the explicit asymptotic behavior of this random matrix. 

In the present paper, we focus on the mono-channel case. Our main result is Theorem 14. 1[ 
It relies on a result obtained by one author in [7J and can be stated as follows: 

Theorem 1.1. Let k be an integer and t be a real number in (0,1). Let be a sequence 
of random channels defined according to Equation ([3]). Then there exists a probability vector 
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(defined in Equation (j4]) j such that, for all e > 0, almost surely when n oo, for all 
input density matrix p, the inequality 

(1) spec(«>(/9)) ^ /JW 

is e-close to being satisfied. Moreover, f3^^^ is optimal in the sense that any other probability 
vector 13 S satisfying the same property must satisfy (3^^^ -< (3. 

We combine this result result with bi-channel bounds to obtain new counterexamples to 
the additivity conjectures for p > 1. 

An illustration of our result is Corollary 15.61 which one can reformulate as follows: 

Theorem 1.2. For each p > 1 and each finite quantum space A of dimension k' ^ 2, there 
exists an integer such that for each quantum system B of dimension larger than this integer, the 
quantum channel arising from a quantum coupling A and B (of appropriate relative dimension, 
depending on p and k' ) has a high probability to be Renyi superadditive when coupled with its 
conjugate. 

From a quantum information theory point of view, the true novelty of this result is that 
any dimension A;' ^ 2 is acceptable. This result does not seem to be attainable with the 
alternative proofs available in [14^ [T2| (Sj [10] . 

Our techniques rely on free probability theory. They allow us to understand entanglement 
of random subspaces, and do not rely on a specific choice of a measure of entanglement. Even 
though the von Neumann entropy is the most natural measure of entanglement in general, 
this subtlety is important as the papers [U [2] imply that all the p ^ 1 Renyi entropies don't 
enclose enough data to fully understand entanglement. 

Our paper is organized as follows. We first recall a few basics and useful results of free 
probability theory of random matrix theoretical flavor in Section [2j In Section [3l we describe 
the random quantum channels we study. Section U] describes the behavior of the eigenvalues 
of the outputs of random channels. In Section O we use results of the previous sections and 
of [8] to obtain new counterexamples to the additivity conjectures. 

2. A REMINDER OF FREE PROBABILITY 

The following is a summary of results contained in [7j , [22] , [23] and [9] . 

2.1. Asymptotic freeness. A non- commutative probability space is an algebra A with unit 
endowed with a tracial state ip. An element of A is called a (non-commutative) random 
variable. In this paper we shall be mostly concerned with the non-commutative probabil- 
ity space of random matrices (A^„(L°°~(r2, P)), E[n~"^ Tr(-)]) (we use the standard notation 

L°°-(fi,P) = np^iLP(fi,P)). 

Let ^1, . . . ,Ak be subalgebras of A having the same unit as A. They are said to be free if 
for all Oj S Aj^ {i = 1, . . . , k) such that (p{ai) = 0, one has 

ip{ai ■■■ak) = 

as soon as ji ^ j2, j2 / Js, . . . ,jk-i 7^ jk- Collections Si, S2, ... of random variables are said 
to be free if the unital subalgebras they generate are free. 

Let (ai,...,afc) be a A:-tuple of selfadjoint random variables and let C{Xi, . . . , Xi^) be 
the free *-algebra of non commutative polynomials on C generated by the k indeterminates 
Xi, . . . , Xf^. The joint distribution of the family {ajj^^-j^ is the linear form 

Ai{ai,...,afe) : C(Xi, . . . ,Xk) C 

P ^ ip{P{ai,. . .,ak)). 
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Given a /c-tuple (ai, . . . , a^) of free random variables such that the distribution of aj is ^q. , 
the joint distribution /i(ai,...,afc) is uniquely determined by the fJ-ai's. A family (a", . . . ,a^)ra 
of /c-tuples of random variables is said to converge in distribution towards (oi, . . . ,0^) iff for 
all P G C(Xi, . . . , Xfc), lJ'[a'^^,„^a]^){P) converges towards M(ai,...,afc)(^) as n — > 00. Sequences 
of random variables (a")n) ■ ■ ■ , iO'k)n called asymptotically free as n ^ 00 iff the /c-tuple 
(a", . . . , a^)n converges in distribution towards a family of free random variables. 

The following result was contained in [22] (see also [9]). 

Theorem 2.1. Let {uj^^}keN be a collection of independent Haar distributed random matrices 
of M.n{C) and {l^^^^^jfcgN be a set of constant matrices of M.n{^) admitting a joint limit 
distribution as n ^ 00 with respect to the state n^-'^Tr. Then the family VF^^^j^g^ 
admits a limit *- distribution {uk,Wk}ke'H with respect to E[n^^Tr], such that ui, U2, 
{wi,W2, ■ ■ ■} are free. 

2.2. Free projectors. Let us fix real numbers ^ a,/5 ^ 1, and consider, for all n, a 
selfadjoint projector 7r„ G A^n(C) of rank qn such that asymptotically qn ~ an as n — > 00. 
Let ir'n be a projector of rank q!^ such that q^^ ~ /3n, and assume that it can be written under 
the form UTTnU* , where C/ is a Haar distributed unitary random matrix independent from iTn 
It is a consequence of Theorem 12. H that -7r„ and vr^ are asymptotically free. Therefore 
TTn'^n'^n has an empirical eigenvalues distribution converging towards a probability measure. 
This measure is usually denoted by //i Kl /i2 , where /ii , //2 are the limit empirical eigenvalue 
distributions of the projectors vr^ and vr^ respectively: 



yUi = (1 — a)6Q + a6i 



= {I - P)5o + P6i. 



In this specific case, we can compute explicitly /xi Kl fj,2- For this purpose, we introduce two 
2-variable functions which will be of great importance in what follows. 

ip+:{{x,y) G [0, ^ [0, 1] 
^ 1- 

if- -.{{x^y) G [0,1]2}^ [0,1 

1 - 

Let us omit the variables of w"^/" and rewrite 



+/- = (/9+/-(a,/3) =a + P-2ap± V4a/3(1 - a)(l - /3) 



It follows then from [23], Example 3.6.7, that 



]dx 



/ii K //2 = [1 - min(a, f3)]6o + [max(a + /? - 1, 0)]6i + 

The proof relies on a technique introduced by Voiculescu to compute ^1 M ^2 in general, called 
the 5'-transform. For more details, we refer the interested reader to [23]. Since we are only 
interested in (^+, we consider the two-variable function ip : [0, 1]^ — > [0, 1] 



V>{x,y) 



if X = or y = 0; 
(p'^{x,y) if X, y > and X + y ^ 1; 

1 if X + y > 1. 
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In the case where a + P < 1, the ranges of vr^ and vr^ do not (genericahy) overlap and 
ip{a,(3) < 1. The previous asymptotic freeness results imply that almost surely, 

liminf ||7r„7r^7r„||oo ^ ^{a,(3). 

n 

We are interested in whether we actually have 

lim||7r„7r^7r„||oo = <f{a,l3) < 1. 

n 

This turns out to be true. This is an involved result whose proof we won't discuss here. We 
just recall the result below as a theorem, following [7] (Theorem 4.15), see also |19) . 

Theorem 2.2. In C", choose at random according to the Haar measure two independent 
subspaces Vn and of respective dimensions qn ~ an and q'^ ~ (3n where a,f3£ (0, 1). Let 
Tin (resp. ir'^) be the orthogonal projection onto Vn (resp. V.^). Then, 

lim||7r„7r.^7r„|| =ip{a,P). 



3. Quantum channels and additivity conjectures 

3.1. Renyi entropies and minimum output entropies. Let = {x S | X]i=i = 
1} be the {k — l)-dimensional probability simplex. For a positive real number p > 0, define 
the Renyi entropy of order p of a probability vector x £ to be 

k 

^ i=l 

Since limp^i H^^x) exists, we define the Shannon entropy of x to be this limit, namely: 

k 

H{x) = H^{x) = - "^XilogXi. 

i=l 

We extend these definitions to density matrices by functional calculus: 

^^(p) = 7^1ogTr/; 
1 — p 

H{p) = H\p) = -i:ip\ogp. 

Given a vector x G C", ||x|| = 1, we call Px the rank one orthogonal projection onto the span 
of X. Using Dirac's bra-ket notation, Px = \x){x\. More generally, for a subspace V C C^, we 
denote by Py the orthogonal projection onto V in A^„(C). 

A quantum channel is a linear completely positive trace preserving map ^ : A4„(C) 
A4fc(C). The trace preservation condition means that density matrices are mapped to density 
matrices, and the complete positivity reads: 

Vd ^ 1, <I> (8) Id : A4„d(C) Mkd{C) is a positive map. 

We recall that according to Stinespring theorem, a linear map ^ : A4„(C) A^jt(C) is a 
quantum channel if and only if there exists a finite dimensional Hilbert space JC = C^, and a 
partial isometry V G End(C"', C''''^) (satisfying V*V = In) such that 



(2) 



^X) = Tv,c[VXV*], yXeMniC). 
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For a quantum channel <I> : A^„(C) A^a;(C), we define its minimum output Renyi entropy 
(of order p) by 

<in(^) = min H^iHp)). 

p£M„{C) 
p^O,Trp=l 

Since the Renyi entropies are concave functions, their minima are attained on the extremal 
points of the set of density matrices and hence 

<in(^) = mm HP{^iP,)). 

\\x\\=l 

3.2. The random quantum channel model. We fix an integer k and a real number t £ 
(0,1). For each n, let Un G Ainki^) be a random unitary matrix distributed according to 
the Haar measure, and (?„ be a projection of 7W„fc(C) of trace pn such that pn/{kn) ~ t 
as n — > oo. To qn we associate a non-unital matrix algebra map Xn '■ -^p„(C) A^„fc(C) 
satisfying Xn(l) = Qn- The choice of Xn is unique up to unitary conjugation, and the actual 
choice of Xn is irrelevant for the computations we want to perform - in the sense that any 
choice will yield the same results. 

We study the sequence of random channels 

given by 

(3) '^n{X)=TTn{Un{Xn{X))U*). 

Remark 3.1. In our previous paper [8j, we considered exactly the same model of random 
quantum channels, with one small difference: the partial trace was taken with respect to C^. 
However, it is well-known that, when partial tracing a rank one projector, the non-zero eigen- 
value of the resulting matrix do not depend on which space is traced out. Hence, from the 
point of view of eigenvalue statistics, the model we consider here is identical with the one in 
[8], Section 6.2. 

Graphically, our model amounts to Figure [H We refer to the first paper of this series, ^ 
for details about this graphical notation. 





c 




k-A 
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-<X' 



Figure 1. Diagram for ^{X) 



We are interested in the random process given by the set of all possible eigenvalues of ^{X) 
as n — > oo. In our setup, we deal with k eigenvalues. 

Let Vn be the image of UxnU*. This is a random vector space in C" (8> of dimension pn 
distributed according to the uniform measure on the Grassmannian space Grp^(C"*^). 

If we can ensure that the entanglement of every norm one vector x G in C"" is 
large with high probability for the uniform measure on G Grp„(C"''^), this will yield new 
entropy bounds. The entanglement is always a concave function of the principal values of x. 
We recall that for an element x S C" we denote A(x) and rk(j;) the singular values and 
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the rank of x, when viewed as a matrix x G MnxkiC). In quantum information theory, these 
quantities are also cahed the Schmidt coefficients and the Schmidt rank of x respectively: 

rk(x) 

X = ^ y^Xi{x)ei (gi fi, 
i=l 

where {cj} and {fi} are orthonormal families from C" and C'^ respectively. Both quantities 
can also be expressed as the rank and respectively the spectrum of the reduced density matrix 
Tr„ Px- The strategy adopted in this paper is to describe a convex polyhedron such that with 
high probability, for a vector subspace V chosen at random, for all input x £ V, the eigenvalue 
vector A(x) belongs to a neighborhood of this convex set. 

3.3. Known bounds. Some results are already available in order to quantify the entan- 
glement of generic spaces in Grp^(C" C'^). The best result known so far is arguably the 
following theorem of Hayden, Leung and Winter in |15j : 

Theorem 3.2 (Hayden, Leung, Winter, [15j . Theorem IV. 1). Let A and B he quantun systems 
of dimesion da and ds with ^ (i^ ^ 3 Let < a < logd^- Then there exists a subspace 
S C A B of dimension 

d ~ dAdsy, T~v)^ 

[logdAy-^ 

such that all states x £ S have entanglement satisfying 

H{P^)^\ogdA-a- 13, 
where (3 = ^^1/(^5 log 2) and T = 1/1753. 

To prove this result, the authors require sophisticated methods from asymptotic geometry 
theory. In particular, they need estimates on the covering numbers of unitary groups by balls 
of radius e and results of concentration of measure. The results of concentration of measure are 
applied to a specific measure of entanglement (e.g. one entropy H^), therefore the measure of 
entanglement does not deal directly with the behavior of the Schmidt coefficients, but rather 
with the behavior of a function of them. 

4. Confining the eigenvalues almost surely 

4.1. Main result. Our strategy is to describe a convex polyhedron K inside the probability 
simplex with the property that, for all e > 0, almost surely when n goes to infinity, all 
input density matrices p ^ 0,Tr(/9) = 1, are mapped to output states <I>(p) whose spectra are 
contained K + e, the e-neighborhood of K. 

For t £ (0, 1), let us first define the vector /3*^*) G M'^', where 

(4) pf) = ^(j.^t^-^(^l_l^ty yi^j^k. 

One can check directly that is a probability vector and that it is a non-increasing sequence. 
Moreover, = (p{l/k,t) and = for j ^ [k{l - t)\ + 2. 

Since the majorization partial order plays an important role in this situation, let us remind 
here the definition and some basic properties of this relation. For two probability vectors 
x,y £ Afc, we say that x is majorized by y (and we write x ^ y) iff for all j £ {1, . . . ,k} 

j j 

(5) Sj (x) = ^ a;} ^ ^ = sj (y) , 

i=l i=l 
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where and are the decreasing rearrangements of x and y; note that for j = k we actually 
have an equality, since x and y are probability vectors. We extend the functions Sj, by 
functional calculus, to selfadjoint matrices X £ A^fc(C). The majorization relation can also 
be characterized in the following way: for a probability vector y and a permutation a G &k, 
denote by a.y the vector obtained by permuting the coordinates of y along a : {a.y)i = y^ij^y 
Then 

a; -< y iff. x £ S{y), 

where S{y) is the convex hull of the set {a.y \ a £ &k}- Moreover, the extremal points of S{y) 
are exactly y and its permutations a.y. In Figure [21 we plot A3, the 2-dimensional simplex 
together with the sets for t = l/k' and k' = 2,3,4,5,10,20,50,100. Notice that for 

k' = 2,3, the set S{I3^^^^ touches the triangle A3, because of the fact that has in this 
case a zero coordinate. 




Figure 2. The 2-dimensional probability simplex with the sets S{P^^^), for 
t = l/k' and k' = 2, 3, 4, 5, 10, 20, 50, 100. 

We can now state the main result of this section: 

Theorem 4.1. Let t be a parameter in (0, 1) and e > 0. Let S'(/3*-*^) + £ be the e-ball around 
S{P^^^) in Afc. Then, almost surely when n 00, for all input density matrix p, 

(6) spec($(p)) G +e. 

Moreover, j3^^^ is optimal: a probability vector [3 G A^ such that, with positive probability, 

(7) spec(^'(p)) G +e Vp 
must satisfy /J^*-* ^ /3. 

We split the proof of Theorem 14. 1 1 into several lemmas. The first one is an easy consequence 
of the definition of the operator norm. 

Lemma 4.2. Let Q,R be two selfadjoint projections in A4„(C). Then 

\\QRQ\\^= meixTriP^R). 

■rczlrn (.) 
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Proof. Since QRQ is a self adjoint operator, we have: 



WQRQW^ = sup {QRQy,y) = sup {RQy,Qy) 

lls/ll^i llyll^i 

= sup {Rx,x) = max Tic{PxR). 

x€lmQ xGlmQ 



The following lemma is a reformulation of the min-max theorem: 



□ 



Lemma 4.3. Let Ai ^ A2 ^ • • • ^ be the Schmidt coefficients of a vector x £ C"^. Then, 
for all 1 ^ j ^ k, 

Sj{x) = Ai + A2 H h A,- = max Ti{PxPc"^f)- 

Proof Since Aj are the eigenvalues of Tr„ Px £ A4k{C), the min-max theorem for Tr„ P^ can 
be stated as: 

Sj{x) = max Ti^PpTrnPx)- 
The conditional expectation property of the partial trace implies that 

Sj{x) = max Ti:{Px ■ In <^Pf) = max Tr(Pa;. • Pc"®^)- 

FGGr,(Cfe) FGGrj(Cfe) 

□ 

We are interested in majorization inequalities which hold uniformly for all norm one ele- 
ments of a subspace V. In other words, we are interested in the quantity 

max Si(x) = max max Tr(PrPrnK(,p). 
x&V ■'^ ' x&V FdGrJC^^) 
\\x\\=l \\x\\=l 

Since A; is a fixed parameter of our model, in order to compute the maximum over the Grass- 
mannian, it suffices to consider only a finite number of subspaces F: 

Lemma 4.4. For all e > 0, for all j, there exists a finite number of j- dimensional subspaces 
Fi,...,Fn € Grj(C*=) such that, for all x £ , 

N N 
m&yiTi{PxP(cr,^Fi) ^ Sj{x) ^ maxTr(P2:Pc"(g)Fj + £■ 

i=l i=l 

Note that in Lemma 14.41 ^ does depend on e but can be chosen to be finite for any e > 0. 

Proof We only need to prove the second inequality. Since the Grassmannian Grj(C'^) is 
compact and metric for d{E,F) = \\Pe — PfWoo^ all e > there exists a covering of 
Grj(C*'') by a finite number of balls of radius e centered in Fi,. . . ,F]\f. Fix some x G C""*^ 
and consider the element F S Grj(C'^) for which the maximum in the definition of Sj{x) is 
attained. F is inside some ball centered at Fi and we have 

Tr(P^.Pc"(g.F) ^ Tr(P^Pc"(g)Fj + |Tr(P^(Pc"(g)F) - Pc"»F,))\ = 

= TriPxPcn^F,) + \\Pf - Pf, IL ^ Tr(P^Pcn^^J + e, 
and the conclusion follows. □ 

Now we are ready to prove Theorem 14.11 
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Proof of Theorem \4-l\ First, notice that it suffices to show ([6]) holds for rank one projectors 
p = Px- The general case follows from the convexity of the functions si, . . . , s^- 

Let e > and j G {1, . . . , k}. For a random subspace V C C"'^ of dimension p„ ~ tnk, 

max sAx) = max max TrfPr-Pr^ssF)- 

||x|| = l ||a::|| = l 

Using the compactness argument in Lemma 14.41 one can consider (at a cost of e) only a finite 
number of subspaces F: 

N 

max Sj{x) ^ max max TT{PxPcn^p.) + e. 

x£V i=l xGV * 

||2;||=1 ||a;||=l 

According to Theorem 12.21 foi" all i G {1, . . . N}, almost surely when n ^ oo, 

lim\\PvPc"^F,Pv\\oo = y^{j/k,t). 

Since N is finite, with probability one, the above equality is true for all i. Next, using Lemma 
14.21 one has that, almost surely, 

limsup max Sj{x) ^ (p{j/k,t) + e, 

n 3:&V 
\\x\\=l 

which concludes the proof of the direct implication. 

Conversely, let fi G A^, be a probability vector which satisfies Equation ([Tj). For j G 
{1,2,..., k} fixed, let Fq be a subspace of C'^' of dimension j. We have 

max sAx) = max max TrfP^. -Pc^ssf) 

x^V ' x&V FeGr,(C*) ^ ^ ' 

||x||=l ||a;||=l 

^ max Tr(P^Pc"®Fo) = WPvPc^^FqPvWoo > ^{j/k,t). 

x£V n-*oo 

Since, with positive probability, max x^v Sj{x) ^ Sj{(3) + e, we conclude that Sj{(3) ^ 
(f{j/k,t) = Sj{P^^^) and the proof is complete. □ 

The interest of Theorem l4.1l in comparison to Theorem l3.2l is that it does not rely specifically 
on one measurement of entanglement, as we are able to confine almost surely the eigenvalues 
in a convex set. Also, our argument relies neither on concentration inequalities nor on net 
estimates, as we fix k. However, unlike Theorem 13.21 Theorem 14.11 does not give explicit 
control on n. It is theoretically possible to give an explicit control on n (using techniques 
introduced in [19]), but this would lead to considerably involved technicalities. 

4.2. Application to Entropies. Once the eigenvalues of the output of a channel have been 
confined inside a fixed convex polyhedron, entropy inequalities follow easily. Indeed, the 
confining polyhedron is defined in terms of the majorization partial order, and thus the notion 
of Schur-convexity (see [1]) is crucial in what follows. 

A function / : M'^ ^ M is said to be Schur-convex if x ^ y implies f{x) ^ f{y)- The Renyi 
entropies H'^ are Schur-concave, and thus majorization relations x < y imply H^{x) ^ iF^{y) 
for all p ^ 1. The reciprocal implication has been studied in [H [2]: entropy inequalities 
HP{x) ^ H^{y) (for all p ^ 1) characterize a weaker form of majorization called catalytic 
majorization, which has applications in LOCC protocols for the transformation of bipartite 
states. 

For the purposes of this paper, the main corollary of Theorem 14. II is the following 
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Theorem 4.5. For a fixed parameter t, almost surely, when n oo, for all input p G T>tnk, 

liminf<iJci,^)^//^'(/3W). 

n 

Proof. This follows directly from Theorem 14.11 and from the Schur-concavity of the Renyi 
entropies. □ 

5. New examples and counterexamples of superadditive channels 

Since our main result, Theorem 14.11 is valid almost surely in the limit n — > oo, the limiting 
objects depend only on the (a priori fixed) parameters k and t. In what follows, we consider 
large values of the parameter and introduce the "little-o" notation o(-) with respect to the 
limit k ^ oo. 

5.1. Superadditivity. We start with a crucial recent series of result, which we summarize 
into the following theorem: 

Theorem 5.1. For all p ^ 1, there exist quantum channels <I>i and ^2 such that 
(8) ® ^2) < Hf^.J'^>^) + if^in($2). 

This theorem results mainly from the papers [121 114^ 116]. Note that the equality 

for any <I>i and ^2 and any p > 1 was a conjecture until 2007, and that even nowadays, no 
concrete counterexamples are known for p small or p = 1. 

5.2. The Bell phenomenon. In order to provide counterexamples for the additivity con- 
jectures, one has to produce lower bounds for the minimum output entropy of single copies 
of the channels (and this is where Theorem 14. II is useful) and upper bounds for the minimum 
output entropy of the tensor product of the quantum channels. The latter task is somewhat 
easier, since one has to exhibit a particular input state such that the output has low entropy. 

The choice of the input state for the product channel is guided by the following observation. 
It is clear that if one chooses a product input state p = pi<^ P2, then the output state is still 
in product form, and the entropies add up: 

HP{[^1 ® ^2]{P1 ® P2)) = HP{^i{pi) $2(P2)) = HP{^,{p,)) + HP{<^2{P2))- 

Hence, such choices cannot violate the additivity of Renyi entropies. Instead, one has to look 
at entangled states, and the maximally entangled states are obvious candidates. 
All our examples rely on the study of the product of conjugate channels 

where 

$„(X) = TlniUnXniX)U:), ¥„(X) = Tr„(I7„Xn(^)f/* ) 

have been introduced in subsection 13.21 Our task is to obtain a good upper bound for 

lim sup i/^. „($„«)¥„). 

n 

Our strategy is systematically to write 

limSUpF^.„($„ ®¥„) < H^^J^n^^niEtnk)) 
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where Etnk is the maximaUy entangled state over the input space (C*"'^)'^^. More precisely, 
Etnk is the projection on the Bell vector 

^ tnk 

Belltnfc = -y= Ci^Ci, 



i=l 



where {cj}*"^ is a fixed basis of C^^^. Using the graphical formalism of [8], we are dealing 
with the diagram in Figure [3] (recall that square symbols correspond to C*', round symbols 
correspond to C", diamond ones to C*"'^ and triangle-shaped symbols correspond to C^^*). 



Z ^ (g> (Etnk) 



u 


t-A ^_ A-t. 

tnk 

» — ^ /♦ — C 


U* 









u 



u* 



u 



u 



u 



u 



Figure 3. Zn = ^'^ ^^{E, 



tnk) 



The random matrix <^n{Etnk) was thoroughly studied in our previous paper [8] and 
we recall here one of the main results of this paper: 

Theorem 5.2. Almost surely, as n ^ oo, the random matrix <I>„(Bellt„fc) G Mi^2{C) 
has eigenvalues 



7 



(t) 



\ 



t + 



1-t I 



1-t 



v 



fc-^ — 1 times / 



From this we deduce the following corollary, which gives an upper bound for the minimum 
output entropy for the product channel $ 



Corollary 5.3. Almost surely, as n ^ oo, 

1 



limsupi?^„($n«'$n) ^ 



■log 



n 1 - p 

In the case p = 1 the upper bound simplifies to 



1 - t\^ 
t + ^] +ik^ 



limSUpFmin(^'n <E)<^n) ^ - {t + 



1 - t 



log U + 



1 - t 



{k- 



1-t 



1-t fl-t 
11^7^ log 



A;2 



A;2 



5.3. Macroscopic counterexamples for the Renyi entropy. In this section, we start 
by fixing t = 1/2. We assume that k is even, in order to avoid non-integer dimensions. A 
value of 1/2 for t means that the environment to which the input of the channel is coupled is 
2-dimensional, i.e. a single qubit. The main result of this section is that we obtain a violation 
of the Renyi entropy in this simplest purely quantum case, k' = 2. 

Using Theorem (521 the asymptotic eigenvalue vector for the output of the product channel 

is 
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The series expansion for H^{'y) when k ^ oo and the Corollary 15.31 imply that, almost surely, 



(9) 



limsupF^i^($ ^>) ^ 



p — 1 



log2 + o(l) 



In the case of a single channel, since ^^ix, 1/2) = 1/2 + y^x(l — x), the vector /? = /?(^/^) is 
has a particularly simple form in this case: 



_i Vfc£T 



i-1 



, V2^j^fc/2, 



Pj = 0, V /i:/2 < i ^ /c. 



where ^'(a;) = x{l — x). Note that the first eigenvalue is large (of order 1/2) and that the 
others are small: 



Theorem 5.4. Almost surely as n ^ oo, 



hminf <;„($)= liminf<ijcl>) ^ 



p — 1 



log2 + o(l). 



Since 



limsupiJ^jjj($ ^>) ^ 



p — 1 



log2 + o(l). 



the additivity of the Renyi p-norms is violated for all p > 1. 

Proof. We shall provide a lower bound for H^{(3). Notice that the main contribution is given 
by the largest eigenvalue: /J^ = 2~^ + o(l). Next, we show that the contribution of the smaller 
eigenvalues is asymptotically zero. We consider three cases: p > 2, p = 2 and 1 < < 2. If 
p > 2, then 

^ j;fc-P/2^A:i-P/2 = o(l). 



For p = 2, one has: 



j=2 
k/2 

i=2 



'^/S r . • _ 1 -I 2 

^(^)-^(V) 



-\ 2 



k/2 

i=2 
fc/2-1 



sup 

0-l)/fc^a;«:i/fc 



tij'ix) 



1 (1 - 2j/kr 

k ^ 4j(l-i/A;) ^ 



The case 1 < p < 2 is more involved: 



k/2 

J ^2 j=2 



1 // 



k 



IP 



1/2 



%b'{t)Pdt 



0(1). 



Hence, in all three cases, HP{[5) ^ log 2 + o(l). This inequality and Eq. ([9]) provide the 
announced violation of the additivity conjecture for Renyi entropies. 



p — 1 



log 2 < 2 • 



P 



p — 1 



log2 + o(l) 



^ 2HP- ($„). 



□ 
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Let us now come back to the more general case of arbitrary t G (0, 1) fixed. It is natural to 
ask whether the bound HP{P^^^) is optimal. Even though this is an open question for fixed k, 
the corollary below implies that it is asymptotically optimal for large k. More precisely, let 
^k,n be the random quantum channel introduced in Section [3.21 (since k will vary in the 
statement below, we need to keep track of it). We can then state the following 

Corollary 5.5. For all p > 1, there exists a sequence tending to infinity as k tends to 
infinity, such that, almost surely 

In particular this means that we can almost surely estimate the Schatten Si — > Sp norm of 
that quantum channel: 

k k 

Proof. For t = 1/2, this follows directly by a diagonal argument from Theorem 15.41 and 
Equation ([9]) together with the simple fact that the entropy increases when one takes tensor 
products: 

The asymptotic estimates of Theorem 15.41 are readily adapted to arbitrary t E (0, 1). As for 
the norm estimate, it follows from the definition of the Schatten norm and the Renyi entropy, 
as well as the fact that the Si Sp norm is attained on density matrices. □ 

It is remarkable that the norm estimate for ||<I> <^||5j^5p given by ||$ (8) ^{Etnk)\\Sp is 
actually optimal. The above corollary stands as a mathematical evidence that the Bell states 
asymptotically maximize the 5*1 — > Sp norm of <I> (8) ^. 

The first example of 'Renyi superadditive' quantum channel was obtained by Holevo and 
Werner in using a deterministic channel. However, their example violated the additivity 
conjecture only for p > 4.79. Hayden and Winter found a class of random counter examples 
for the whole range of parameters p > 1 m |14j . Our being able to prescribe t in the coun- 
terexample of Theorem 15.41 is an improvement to the counterexamples provided in the paper 
[T6] (even though there is evidence that the very recent techniques of [IHl 15] could be applied 
for p > 1 and finite t - yet perhaps not as big as 1/2 or 1/3). 

Physically, this means that to obtain a counterexample, it is enough to couple randomly 
the input to a qubit (A;' = 1/t = 2) to obtain a counterexample. The above reasoning applies 
actually for any t. In the following corollary we focus on the case t = 1/k' for integer k', as 
it is more relevant physically. 

Corollary 5.6. For each p > 1 and each integer k' ^ 2, let t = 1/k'. There exists an integer 
kf) such that for all k ^ kQ, one has almost surely 

limsupi?^^($(g)¥) < 21iminf i7^^(^>) 

n " 

Since the proof is very similar to the case k' = 2, instead of providing the details, we plot 
in Figure m acceptable values for feg as functions of p, for several values of A;' = 1/t: 

ko{t,p) = mm{k G N | lim sup HP{<^ (^^){Etnk) < 2FP(/3W)}. 

n 

Note that fco as defined above may note be the smallest dimension yielding a violation of 
p-Renyi additivity. It may be that a better choice for the input state of the product channel 
could yield a smaller value for H^^^{^ $). As the plots suggest, the values of k^ are not 
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bounded when p — > 1. This fact is independent on the choice of the parameter t = 1/k' . The 
results of [5l [10] suggest that there should be a k' large enough for which it is possible to keep 
ko bounded as p — > 1. This improvement is due to their better bounds on H^^j^{^), obtained 
using the techniques developed by Hastings in [12] . 




24 



22 



20 



k_0 



18 



16- 



„1 1,2 1,3 1,4 1,5 1,6 1,7 1,8 1,9 2,0 

p 



(a) 



(b) 



k'=2 
-k'-3 
■ k'-lO 



Figure 4. Plots of ko, i.e. acceptable values of k for which we get an asymp- 
totic additivity violation, in function of p, for different values of k' = 1/t. Two 
ranges for p are plotted separately: p £ [1.1, 2] in (a) and p £ [2, 5] in (b). 



We finish this section by a computation showing that the above bounds are not good 
enough to obtain the violation of the additivity conjecture in the case p = 1. We start with 
the entropy of the product channel: 

F(7) = logfc + log2 + o(l). 

For the case of the single channel, we need an upper bound for H{P^^^'^^) (recall that h{x) = 
—xlogx): 

k/2 k/2 



i=2 

k/2 

i=2 



i=2 



log 2 



sup ^'{x) 



{j-l)/k^x^/k 



+ o{l) 



log 2 



+ o{l) 

k/2 



i=2 



1 

■0 



A;" V k 



Using 



and 



log 2 

[ho^']{t) dt = 



+ 0(1). 
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we obtain 



log2 , logk ( J 



i^(/j'"=')<T+ET*'(i)+Si[''»v/i(f)+o(i 



log 2 1 log 2 log A; 

= ^ + 2 log A: ^ + o(l) = + o(l). 

At the end, the entropy deficit is 

H{j)-2H{I3) ^ log2 + o(l) > 0, 
which does not yield a violation of the minimum output von Neumann entropy. 

5.4. The case t = k~°^. We conclude this paper with the study of the case t = k~" , where 
Q > is a fixed parameter. This corresponds to an exploration of a larger environment 
size C*^". To simplify the computations, we consider only the case of the minimum output 
von Neumann entropy. As before, we provide estimates, when k is fixed but large, for the 
minimum output entropies of $ ® $ and 

We start with the simpler case of the product channel <I> <I>. Theorem 15.21 provides the 
almost sure eigenvalues of ® <I>](£'t„fc): 



7 = 7 



1 _^ 1 



1 1 



1 



1 



1 



v 



Ua ' U2 Ua+2 ' 1.2 h.a+2 ' " " " ' 1,2 h.a+2 

i\j i\j r\j i\j i\j i\j i\j 
^ V 

— 1 times 



Using the series expansion h{l — x) = x — x^/2 + o{x'^), one can compute the asymptotics for 
the minimum output entropy: 

Proposition 5.7. For the product channel <I> Cgi <I>, the following upper bounds hold almost 
surely: 



(2logk-^2^^ + o('^) ^/ 0<a<2; 



(10) ifmin($ $) ^ F(7) = <^ 



21og/c 



2 log 2-1 



if a = 2; 



2\ogk- 2^'L-'i + o (p^) if a>2. 
Our estimate for the single channel case is as follows: 
Proposition 5.8. For all a > 0, the following lower bound holds true almost surely: 



log k 



+ 



log k 



For the purposes of this proof, we define : [0, 1 — /c "] ^ [0, 1], ^k{x) = ^{x^ k ") and 
h{x) = — xlogx. We have 



d_ 

dx 



fix, y) = l-2y + yjy{l - y) 



where the function g : (0, 1) 



X - X 

is defined by 

\/l — X y/x 



(1 - 2?/) 1 + 



gjx) 
9{y) 



g{x) 
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Proof. According to Theorem 14.51 for e > 0, 

where P = P^'' is the /c-dimensional vector defined by 



/3i 



1 / (ij 



V 2 ^ j ^ J; 



y J <j f^k. 



The index J is the number of non-trivial inequalities we get by using Theorem 14.11 and it 
is equal to A: — 1 if a ^ 1 and to [k — /c^~"J if a < 1. 

Our purpose in what follows is to provide a "good" estimate for H{P). We start by rescaling 
the eigenvalues: H[P) = log A: + ^H{kp). In this way, we can focus on the "entropy defect" 
\ogk — H{P) and reduce our problem to showing that 



k 



a-l 



log k 



H{kp) 



k 



log A; ^ 

j=i 



h{kpj 



-1, 



The next step in our asymptotic computation is to replace the unknown points by 
simpler estimates of the type j/k. Notice that the largest eigenvalue /3i is of order k~^. By 
the continuity of the function h, there exists a constant C > such that \h{kPj)\ ^ C and 
thus, individual terms in the sum (llip have no asymptotic contribution. Moreover, we can 
assume J = k — 1, ignoring at most k^~°' terms which have again no asymptotic contribution. 
It is clear that the function x i— > f'i^{x) is decreasing at fixed k and since the entropy function 
h is increasing for x G (0, e~^) and decreasing for x ^ e~^,we can bound h{ip'j^{^j / k)) by 
h{ip'i^{j / k)) , and we reduce our problem to showing that 



(12) 

or, equivalently. 



.a~l 



fe-1 



log A; ^ 



^1. 



-1, 



log A; 



Now, 



h 



(1 - 22/) 1 + 



The term 2h{l - 
h{l + t) + h{l-t) 



9{y) 

= 2h{l - 2y) + (1 - 2y) 
2y) = 2h{l - 2A:-") - 



+ h 



h[l + 



:i-2y) 1 



9{x) 

9{y) 



+ h{ 1 



g{x) 
9{y) 

g{x 



-1. 



9{y) 



~ 4A; " has no asymptotic contribution and, using 
-t^ + 0{t'^), we are left with computing the limit of the main contribution 



k^ 



log k 



1 Lfc/2J 

E 



9ij/kf 

9{yY 



k 



lk/2\ 

y 

log k ^-^ 



-9{jlkfk-\ 
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Finally, 

A; log A; ^ ^ ' ' logA; ^ j l-j/k log A; fc^oo' 
The error term 

, „_i Lfc/2J . . .4 , Lfc/2J ^ Lfc/2J ^ ,^ „ . .4 

V- 9{j/kr fc" _„ri7A-^4,-2a 1 V- 1 (1 - 2j/fc)^ 

logA; ^ logA; ^ ^^-^^ ^ fc^+ilogA; ^ j2 (l-jyA;)2 

converges to zero. In conclusion, we have shown that Equation ([IT]) holds and we deduce that 

rjfn^ 1 , log A; flogk 
H{(3) = log A; - — - + ' 



A;° V , 

□ 

The bounds obtained in this section do not yield a violation of the Holevo additivity con- 
jecture. However, after the first version of this paper was released, Brandao-Horodecki [5] 
and Pukuda-King [TO] used the same model as ours and adapted original ideas from Hastings 
|12j to prove that this model can also lead to a violation of the minimum output entropy 
additivity. 

The techniques in [5l [10] yield more information on the possibility of large values of the 
minimum output entropy for the model under discussion. However, our proofs are of free 
probabilistic nature and yield results of almost sure nature. In addition, [5l [10] rely very 
much on the actual properties of Shannon's entropy function, whereas our techniques attack 
directly the question of the behavior of the eigenvalues. 

We conjecture that the set S{P^^^) (having the property that for any e > 0, S{(3^^^) + e 
contains almost surely the eigenvalues of outputs of random quantum channels) can be made 
smaller and actually optimal, thus yielding as a byproduct that all the values -ffmin(^) converge 
almost surely. However, the results of this paper show that the notion of majorization is not 
sufficient to achieve this goal. 
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